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a transition metal element. We focus on the competition between Coulomb interaction, spin-orbit 
coupling, and lattice distortion when these energy scales are comparable. Strong spin-orbit coupling 
entangles the spin and the t2g d-orbitals giving rise to doublet j — 1/2 and quadruplet j = 3/2 states. 
In contrast to previous works which focused on the doublet manifold, we also discuss the quadruplet 
manifold which is relevant for several pyrochlore oxides. The Coulomb interaction is taken into 
account by use of the slave-rotor mean field theory and different classes of lattice distortions which 
further split the levels of the quadruplet j = 3/2 manifold are studied. Various topological phases 
are predicted, including exotic strong and weak topological Mott insulating phases. We discuss the 
general structure of the phase diagram for several values of d-shell filling and various symmetry 
classes of lattice distortions. Our results are relevant to the search for exotic topological insulators 
and quantum spin liquids in strongly correlated materials with strong spin-orbit coupling. 
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I. INTRODUCTION 

Transition metal oxid es h ave been an active topic 
of research for decades.!^ In particular, the inter- 
play between Coulomb interaction, spin-orbit coupling, 
and lattice degrees of freedom have made transition 
metal oxides an ideal playground to test new theo- 
ries and discover exotic behaviors. Notable among 
them are high temperature superconductivity,'^ colos- 
sal magnetoresistance,'^ and heavy fermion physics,f^'t2l 
among many other possiblities. ^^'^' The rather localized 
nature of 3d orbitals (compared to M and 5d) in some 
transition metal oxides enhances the on-site electron- 
electron int erac tion and typically makes it a dominant 
energy scale.l^In the 3c? transition metal oxides other in- 
teractions such as spin-orbit and electron-lattice coupling 
are typically small compared to the on-site Coulomb in- 
teraction and ground states with antiferromagnetic order 
are typicaL^iSl 

However, Ad and 5d orbitals in layered perovskites such 
as Sr2Ru04, Sr2Rh04, Sr2lr04, Na2lr03, and the hyper- 
kagome Na4lr308 are more spatially extended and thus 
the Coulomb interaction is typically weaker than those 
with 3d orbitals.EI The more extended nature of the 4c? 
and 5c? orbitals compared to the 3d orbitals leads to a 
greater level splitting in a crystal field, and enhances their 
sensitivity to lattice distortions. In many oxides, the 
transition ions are surrounded by an octahedron of oxy- 
gen atoms, MOg, where M represents a transition metal 
ion. The crystal field splits the 5 degenerate (neglecting 
spin for the moment) d-orbitals into two manifolds (see 
Figjlj;): a lower lying {dxy,dyz,d2x) manifold and a 
higher lying Cg [d^z^-r'^ ,dx^-y^) manifold.^ ^ The energy 
separation between the t2g and Cg levels is conventionally 
denoted "lODq" and is typically on the order of ^^1-4 eV, 
which is large compared to many 3c? compounds .'i^l 



Besides the crystal field, the relativistic spin-orbit cou- 
pling is another energy scale that results from the large 
atomic numbers of heavy transition elements. While in 
the absence of spin-orbit coupling the on-site Coulom b 
interaction is of the same order as the band width,'i^'Sl in- 
clusion of strong spin-orbit coupling modifies the relative 
energy scales Thus, for materials with 4c? and partic- 
ularly 5c? electrons, one expects the appearance of novel 
phases with unconventional electronic structure due to 
the characteristic energy of spin-orbit coupling approach- 
ing that of the Coulomb interactions.'^ 

In a cubic environment, the L = 2 orbital angular mo- 
mentum of the c?-orbitals is projected down to an effec- 
tive angular momentum I = 1 (with a minus sign) in 
the t2g manifold.E^ When the spin-orbit coupling is also 
strong, neither spin nor orbital angular momentum is a 
good quantum number. Instead, the total angular mo- 
mentum, i.e. j = I + s, is a conserved quantity, where 
s is the spin of the electron. Thus, spin-orbit coupling 
splits the t2g orbitals with spin into a j = 1/2 doublet 
and j = 3/2 quadruplet separated by an energy gap pro- 
portional to the strength of the spin-orbit coupling, A. 
(See Figjlj;.) For large enough spin-orbit coupling, the 
new effective spin states lead to a great deal of novel 
Mott insulating states^^^'^I] possible spin liquids in the 
hyper-kagome lattice,'22H29l orbital-oriented exchange cou- 
pling in Kitaev-type models,'^^'^ Dirac semi-metal with 
Fermi arcs,'^'^ the quantum spin Hall effect, '^^ topological 
Mott insulators ,1^^ topological magnetic insulators with 
axionic excitations,'22l25| and possibly high temperature 
superconductivity.'SSl 

Of particular interest in this paper are the topolog- 
ical phases that occur in weak to moderately strongly 
interacting systems with strong spin-orbit coupling. "^^ -^ 
(For excellent recent reviews see Refs. [52|53j .) The pos- 
sible time-reversal invariant topological phases of matter 
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have non-trivial topological features in their global band 
structure and robust edge (surface) states. When strong 
electron correlations are taken into account, spin-orbit 
coupling can give rise to a topological Mott insulating 
phase in which the charge degrees of freedom are com- 
pletely gapped (even on the surface) , but where the spin 
degrees of freedom inherit the non-trivial band topology 
from the weakly interacting limit. ^'"^ Thus, the spin 
degrees of freedom form gapless spin-only edge (surface) 
modes. Such edge (surface) states can in principle be 
detected in thermal transport measurements (but not so 
readily in spin transport as the spin-orbit coupling gener- 
ically destroys all spin conservation laws). 

In this work we focus on the interplay and competi- 
tion between strong correlation effects, spin-orbit cou- 
pling, and lattice distortion that is expected to be im- 
portant in heavy transition metal oxides. In the heavy 
transition metal oxides one expects both the spin-orbit 
couplin g'^ •^^^ and the lattice distortion energies^ to be 
of the order of 0.05 — 0.5 eV, while the interaction energy 
is typically at the higher end of this scale to somewhat 
larger, 0.5 — 2 eV.'^^ "^'^ While the phase diagram of an 
interacting undistorted pyrochlore model with j = 1/2 
has already been studied,!^ we expand those results to 
include the effects of distortions of the local octahedra 
on the phase diagram. We also investigate pyrochlore 
oxides at different d-level fillings with the Fermi energy 
lying in the quadruplet j — 5/2 manifold, which has not 
been considered in previous works. One of our motiva- 
tions is to see if the j = 3/2 manifold can also realize 
the interesting Mott phases of the j — 1/2 manifold.^- 
We find that, indeed, these exotic phases can be realized 
for the j — 3/2 manifold. Moreover, we find that for 
the J = 1/2 manifold "weak" topological variants of the 
exotic Mott phases can also appear in the phase diagram 
when certain types of lattice distortion are present. 

This paper is organized as follows. In Sec. |ll]we derive 
an effective nearest-neighbor tight-binding Hamiltonian 
that properly captures the non-interacting limit of the 
physics, and include an on-site Hubbard interaction term 
to describe electron correlations. In the absence of inter- 
actions and distortions, the ground state is metallic for 
weak spin-orbit coupling, and becomes a strong topolog- 
ical insulator as spin-orbit coupling grows. In Sees. Ill 
|IV[ and |V] we study the effects of Coulomb interaction 
and lattice distortion o n equ al footing using the slave- 
rotor mean field theory.l^SEll As the strong correlation 
limit is approached, a Mott transition occurs and exotic 
phases are realized. We discuss the conditions that favor 
these unusual phases. Finally, we conclude in Sec |VII| 
and outline interesting topics for further study. 



II. DERIVATION OF THE EFFECTIVE 
HAMILTONIAN 

In this work, we restrict our attention to pyrochlore 
oxides of the form A2M2O7, where A is a rare earth el- 
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FIG. 1. (Color online) (a) An illustration of the pyrochlore 
lattice which is composed of corner sharing tetrahedra. Tran- 
sition elements are indicated by black solid circles, (b) Each 
transition ion is surrounded by an oxygen octahedron shown 
by six solid blue (dark grey) circles. A transition ion is located 
at the origin of the local coordinate and is shown in black. We 
study a trigonal distortion preserving C3 symmetry applied 
along the [111] direction (or its equivalent), shown by two yel- 
low (grey) faces, and an elongation preserving C4 symmetry 
along the z-axis of the local coordinate, (c) A schematic rep- 
resentation of the splitting of the bare atomic d-levels (1), due 
to a cubic crystal field arising from the octahedral environ- 
ment (2), unquenched spin-orbit coupling in the t2g manifold 
(3), and a distortion of the octahedron (4). The values of the 
splittings in (4) depend on A and A3_4. 



ement and M is a transition metal element. Examples 
include A2lr207 (A=Y, Pr, Eu or other rear-earth ele- 
ments), Cd20s207, and Cd2Re207. In these materials, 
the transition metal elements form a pyrochlore lattice 
and each M sits in the center of an oxygen octahedron.l^ 
The relevant geometry and coordinate system we use, 
along with the important level splittings, are shown in 
Fig.[l] 

For strong spin-orbit coupling and with cPi^a < 4, 
i.e., less than 4 electrons) the upper j = 1/2 manifold 
is empty and the important electronic structure is given 
by the lower j = 3/2 manifold. In ordered double per- 
ovskites with a local quadruplet, strong on-site interac- 
tions add bi-quadratic and bi-cubic exchange interactions 
to the effective spin exchange Hamiltonian deep in the 
Mott regime.®^ It is also argued that the same local C3 
distortion that we consider in subsequent sections maps 
the exchange Hamiltonian into a pseduo-spin-1/2 model 
that favors magnetic order.l^ 

Here we focus on the w eak to moderately strongly in- 
teraction regimes .'^^'^^'^ We begin by deriving an effec- 
tive tight-binding Hamiltonian between transition ions 
located at the vertices of the corner-shared tetrahedrons 
in the pyrochlore lattice shown in Figjl^. The transi- 
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tion ions we consider have 5d orbitals with three or five 
electrons in the triply degenerate (neglecting spin) t2g 
manifold. The spin-orbit coupling in this manifold has 
the following forrrfS 

Hso = -XI ■ s, (1) 

where I — 1 and s = 1/2 describe the orbital and spin 
degrees of freedom, and A > parameterizes the strength 
of the spin-orbit coupling. That the t2g orbitals can be 
effectively described by angular momentum / = 1 comes 
from the projection of the d-orbit al ang ular momentum 
into the local basis of t2g manifold. '^SEB 

To study the effects of lattice deformations)^ we as- 
sume that the octahedron surrounding an ion can be dis- 
torted in two ways: (1) a trigonal distortion preserving 
local C3 symmetry and (2) an elongation (expansion) of 
octahedra preserving local C4 symmetry. (See FigjljD.) 
The former has been argued to be rather common and 
can be described by the following Hamiltonian on each 
transition metal ion site!^ 

Htri = -A^idl^Jzx + dl^d^y + dl^dxy) + h.c, (2) 

where A3 parameterizes the strength and sign of the C3 
preserving distortion, and the C4 elongation/contraction 
splitting is described bj^^ 

Hei = = A4{nyz + n^j;), (3) 

where A4 parametrizes the strength and sign of the dis- 
tortion, and Iz is the z component of the effective angular 
momentum of the t2g orbitals related to the occupation 
of the dxy orbital by Uxy = Ud— (Iz)^ which follows from 
the constraint Ud = Uxy+nyz+Uzx- {See also Appendix 
[A| ) For an elongation of the tetrahedron, A4 < 0, and 
for a compression of the tetrahedron, A4 > 0. Trigonal 
distortions appear to be more common in real materi- 
als, and the magnitude of the energy splittings can be 
crudely estimated from density functional theory calcu- 
lations based on X-ray determined positions of oxygen 
atoms around the transition metals. We are not aware 
of detailed calculations of this type for the Ad and 5d py- 
rochlore oxides, but closely related 3d systems appear to 
have splittings on the level of 0.01 — 0.5 eV.I^ We take 
this as crude estimate, with the larger end of the energy 
scale probably more likely for the more extended Ad and 
5d orbitals. 

Thus, the local Hamiltonian describing the t2g orbitals 
on each site is 

Hlocal — Hso + Htri + H^l ■ (4) 

The Hamiltonian Q can be easily diagonalized and its 
eigenvectors describe a projection onto the spin-orbit 
plus distortion basis. We will denote the projection by 
a matrix Af , which contains all the information about 
the spin-orbit coupling and the distortion of the octa- 
hedra (all assumed identical so translational invariance 
is preserved). Moreover, due to the presence of time- 
reversal symmetry, the eigenvectors form Kramers pairs. 



A schematic representation of splitting t2g upon includ- 
ing the terms in Eq.Q is shown in Fig{l|:. 

We now turn to a derivation of the effective Hamilto- 
nian. We first assume A = A3 — A4 — 0, i.e. neglect 
the contributions in Q. To obtain the kinetic terms of 
the Hamiltonian, we need to describe the t2g orbitals of 
a single ion in the local coordinate system defined by the 
octahedron of oxygen atoms surrounded the ion, and we 
need the p-orbitals of oxygen in the global coordinate sys- 
tem. The hopping of electrons from one transition metal 
ion to a nearest-neighbor transition metal ion is mediated 
by the oxygen p-orbitals. (We note that for the relatively 
extended 5(i orbitals direct overlap may also be impor- 
tant, as well as further neighbor hopping. ^°) We thus 
compute the p-d overlaps to determine the hopping ma- 
trix elements. The local and g lobal axes are related by 
a set of rotation matrices.'^^^ The combination of rota- 
tion matrices and d-p overlaps gives rise to the following 
Hamiltonian: 

Hd = £d d\^^di~,a + t T^^ Ytr'^l'jadi'Ycr' ' 

(5) 

where i, 7, and a in the sums run over lattice sites, t2g 
orbitals {xy, yz, zx), and spin degrees of freedom, respec- 
tively. The Ed stands for the onsite energy of the degen- 

erate t2g orbitals, and t = ^ is the unrotated hopping 
amplitude depending on the overlap integral Ypd-x and the 
energy difference between p and d orbitals. The parame- 
ter t sets the basic hopping energy scale in the problem. 
Without loss of generality we set Sd — 0. 

The effect of spin-orbit coupling and distortion are in- 
cluded via the projection of Hamiltonian in EqJ5| into 
the eigenvectors of the local Hamiltonian in Eq. (WTusing 
matrix M as follows: 

Ho =J2^c,cl^^^a+t J2 r^\'a'cLc«'a', (6) 
ia <ia,i' a' > 

where Va {a = 1,...,6) stands for the six eigenvalues 
of local Hamiltonian Q, and the matrix F describes 
the hopping between sites given in the local basis via 
F = M*TM^. The 4^(cia) is the creation (annihila- 
tion) operator of an electron at site i and in local state 
a. Finally, we add a Coulomb interaction to obtain 

H = Ho + ^Y.(T.4o.c^.-nd)^ (7) 

i a 

where U is the on-site Coulomb interaction and Ud is the 
number of electrons on the 5d orbital of the transition 
metal ion. In the remainder of this paper, we investigate 
the zero-temperature phase diagram of the full Hamil- 
tonian ([7]), which includes the spin-orbit coupling and 
lattice distortions in Q. 

Before closing this section, it is instructive to take a 
look at the non-interacting limit of our model. In the 
absence of distortion, small spin-orbit coupling favors a 
metallic state for all fillings we consider (because they 
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correspond to partially filled bands). However, strong 
spin-orbit coupling opens a gap at \/t = 2.8 for filling 
rid = and for X/t = 2.5 at filling ti^ = 3 (see FigjsJ. 
The two fillings correspond to the j — 1/2 and j = 3/2 
manifolds in the bands, with one hole/ion in each case. 
Therefore, strong spin-orbit coupling can turn the metal- 
lic band structure arising from the quadruplet j = 3/2 
manifold into a band insulating phase. 

By using the Fu and Kane*'^ construction, we can de- 
termine whether the insulating phase is trivial or topo- 
logical. We make use of the inversion symmetry of the 
model and look for the parity of eigenstates at the time 
reversal invariant momenta. Those parities are related 
to the strong topological Z2 index via the following 
relation, 

(-ir=nn^2™(r.), (8) 

m i 

where the first product is taken over filled bands, the sec- 
ond one over the time reversal invariant momenta F^, and 
Cm(ri) is the corresponding parity eigenvalue of band m 
at time-reversal invariant momentum F^. (Note that only 
one band from each set of Kramer's pairs is included in 
the product. Note also, that we have followed convention 
and used Fj with a single subscript for the time-reversal 
invariant momenta. It should not be confused with the 
multi- indexed F" ^, in ^ that describes the hopping.) 
The index i is a collective index, i.e. i € {ii,i2:*3}: in 
terms of reciprocal lattice vectors K of the pyrochlore lat- 
tice: F,j = l/2(iiKi-fz2K2 -1-13X3) so that 11,12-,^ = 0,1. 
The week indices are defined in a similar way to the 
strong index: 

(-i)-^^-=n n ^2™(r.), (9) 

where j = 1,2,3. Thus, the index (i^o; i^ii^2'^3) defines 
sixteen classes of band insulators. If fo = 1, the state 
is said to be a strong topological insulator (STI) and it 
has time-reversal symmetry protected gapless boundary 
excitations described by an odd number of Dirac cones 
in the surface state Brillouin zone.l^ On the other hand, 
if t'o = but Vj ^ for at least one j £ (1, 2, 3) then the 
state is said to be a weak topological insulator (WTI). 
In this case, gapless surface modes may be present in a 
clean system, but can be destroyed with disorder. As 
we will see, we find both STI and WTI (including strong 



correlation generalizations) in our model when lattice dis- 
tortions are present. (See Fig[2] lower panel.) 

III. SLAVE-ROTOR MEAN FIELD THEORY 
AND SELF-CONSISTENT EQUATIONS 

In this section we apply the slave-rot or me an-field the- 
ory developed by Florens and George^^^EH to treat the 
effect of weak to intermediate strength Coulomb interac- 
tions in the regime where the charge fluctuations remain 
important. In this theory each electron operator is repre- 
sented in terms of a collective phase, conjugate to charge, 
called a rotor and an auxiliary fermion called a spinon as 

where Cia is the electron destruction operator at site i 
with quantum number a, representing the states in ([7|. 
The factor e*^' acting on the charge sector is a rotor low- 
ering operator (with Oi a bosonic field), and fia is the 
fermionic spinon operator. The product of the two re- 
sults in an object with fermi statistics, need for the elec- 
tron. Note the rotor part only carries the charge degree 
of freedom while the spinon part carries the remaining 
degrees of freedom a. Therefore, an electron has natural 
spin-charge separation if a is spin in this representation. 

A constraint should be imposed to retain the physical 
states as 

U + Y^flhc^na, (11) 

where L,; is number of rotors. Using this representation, 
the interacting Hamiltonian ([7| can be written as, 

+ |E^^ (12) 

i 

where a chemical potential fi has been introduced. In 
order to treat the phase 9 and angular momentum L on 
an equal footing, we need to switch from (9, L) to fields 
(6',i9e)5a so that idr9 = ff , which gives L = {i/U)dr9. 
The corresponding action is 

dT[~iLdr9 + H + pdrf]. (13) 

Jo 

We next introduce bosonic operators Xi = e*^' and recast 
the action into 



dT 



/L(5r+^^a-M-/^^)/.a+5][^a,x;a,x,+^(x;a,x,-c.c)+p,(|x,|2-l)]+^ 



<i,i'> 



rVl \r 



-^i' fiafi'c 

(14) 
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where hi and pi are Lagrange multipliers imposing the 
constraints Li + J^a fLfia = and jX^p = 1, respec- 
tively on each site. We have effectively carried out the 
integration over L by using the relation L — {i/U)drd- 
The action ( 14 1 describes the coupled spinon and rotor 
degrees of freedom. We will assume that translational 



symmetry is preserved and decompose ( 14 ) into two parts 



by use of the following mean-field ansatz: 



Qe 



Qf^{x:x,,), 

y ] ^a,a' ifiafi'a') 



(15) 



Then, the action in Eq.( 14 1 can be written a,s S = Sf+Se, 
in which 



dr 



<i,i'> 



(16) 



(17) 



The chief benefit of the above actions Sf and Se is that 
they are quadratic in spinon and rotor fields, and there- 
fore the calculation of the corresponding Green's function 
is straightforward. One can simply use Fourier transfor- 
mation and go to the eigenfunction basis to obtain 
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knj 



{iujn + ej{k))fk„j] 



(18) 



k,n,j 



where k is the momentum, j labels the four sites in a 
unit cell as well as the effective spin degrees of freedom, 
and the dispersion of band j is given by ej{k). The iw„ 

are fermionic Matsubara frequencies and fknj is a linear 
combination of fka that diagonalizes the spinon part of 
the Hamiltonian: 



<i,i'> 



(19) 



The rotor action reads 



X*{k,i^n) 



-^+p + tQeak) 



Xj{k,Vn) 



(20) 



where the parameter U has been replaced by U /2, so 
that the atomic limit is treated correctly.l^ Note that 
we have s et h = 0, which guaranties that the con- 
straint Eq.( 11 1 is satisfied on the mean- field level.I^The 



Vn are bosonic Matsubara frequencies and ^(fc) is re- 
lated to the spectrum of the rotor HamiltonianlS via 

and 



Ci,2(fc) = 2(^l±,^l + ^^<^cos(V2)cos(fc,/2)^ 

^3 4 = —2. In the sum under the radical, p, v run over 
three values and give max{|^i|}=6. 



r 



The spinon and rotor Green's function can be readily 
written as 



Gf{k,iujn) 
Geik.iv) = -2- 



1 



iwn + e(fc) ' 
1 



P + tQei{k) 



(21) 
(22) 



which can be used to to determine the character of spinon 
and rotor excitations. We note that the spinous are 
single-particle-like excitations, while the rotors represent 
collective excitations. 

The self-consistent equations that should be solved to 
determine the phase boundaries in the phase diagram are 
(where the rotors are assumed condensed) 



1 = \xy = 



1 1 

4iV~ 



E 



1 



P + tQeij{k) 



, (23) 



1 1 



24iV/3 



(24) 



E 

n,3,k IT +P + tQe^j{k) 



^nd = j^J2®(^'-'om, (26) 



where N is the number of unit cells in the lattice, /3 is 
the inverse temperature, and @(x) is the step function. 
By use of the Matsubara sunP^ 1/ p'^^G0{k,ii'n) — 
— , ^ the first two equations above in the set 
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starting with ( 23 ) can be simplified a^^ 



1 = 



47V 



u 



Qf = - 



2AN 



4{p 



U 



(27) 



(28) 



The rotor condensed phase is characterized by a 
nonzero value oi Z = (e*^), where the electron opera- 
tor is proportional to the spinon operator: Cia = Zfi^- 
The condition for conden sation is that the gap of the ro- 
tor's spectrum Ag — 2yJU{p + 6tQg) closes. Therefore, 
right at the phase boundary p — —QtQg. Combined with 
equations Eq. ( 27 1 and Eq. ( [28| , Q j and the critical U can 
be determined! ^° ' °"^ l 



IV. 



j = 1/2-BAND MODEL: EFFECT OF 
TRIGONAL DISTORTION 



In this section we discuss the slave-rotor mean-field 
phase diagram of the j = 1/2 model studied in Ref.[T5] 
when a C3 symmetry preserving trigonal distortion of the 
oxygen octahedra is included.^" The effect of the trigo- 
nal distortion on the atomic levels is given by Eq. ([2]) , 
which describes the compression or expansion of the oc- 
tahedra along the [111] direction or any equivalent direc- 
tion in the local coordinate. Fig. [T]d shows a schematic 
of this deformation as indicated by the shaded faces of 
the octagon. This distortion splits the t2g manifold into 
a singlet \a) with energy = —2 A3 and doublet je^) 
with energy Eg — A3. 

The effect of spin-orbit coupling on the doublet can 
be understood by noting that the spin-orbit coupling in 
the subspace spanned by the doublet states acts like a 
Zeeman fieldp^ 



A 



(29) 



where and a act on the pseudospin space spanned by 
e'g states and real spin, respectively. The unit vector n 
points in the direction of the trigonal distortion. Thus a 
gap may open by tuning the spin-orbit coupling in the 
presence of a strong distortion, which is consistent with 
density functional calculations for the iridate Na2lr03.^° 
As discussed in Ref.|60j, in the non-interacting limit a 
strong trigonal distortion can turn a strong topological 
insulator into a metal. We would like to understand to 
what degree this happens in the presence of interactions, 
and to what degree the undistorted, interacting phase di- 
agram of Pesin and Balent^^ is changed by distortions. 
We thus study both interaction and distortion on equal 
footing. After verifying that our calculations successfully 
reproduce the phase diagram in Ref.|T5], we first con- 
sider the case of A3 > 0. Our results are shown in Figj2] 
The upper panel exhibits the phase diagram with posi- 
tive distortion. The thick dashed line separates the rotor 



condensed phase (below the line) from the uncondensed 
part (above the line). At the non-interacting level and 
for weak spin-orbit coupling, the system is in the metallic 
phase. However, a small window of gap opening exists for 
A « 2.8t — 3.3i, which is a strong topological insulating 
(STI) phase with Z2 invariant (1;000). This small win- 
dow forms a narrow gapped region in the phase diagram 
along the J7 = axis. Note the "re-entrant" metallic 
phase in the presence of distortion is different from the 
robust insulating phase found in Ref.[T5] for large A, and 
is qualitatively similar to the non-interacting distortion 
results found in Ref. [50] . 

Within the slave-rotor mean-field theory, one finds that 
a narrow window of STI persists to interactions of order 
the bandwidth, after which it becomes a tiny sliver. A 
metallic phase has mostly replaced what would be the 
STI in the absence of distortions. As interactions are 
further increased, and a Mott transition occurs to gap 
out the rotor degree of freedom (above the dashed line), 
one finds the metal is converted into a gapless Mott in- 
sulator (GMI), which is a type of spin liquid with gapless 
bulk spin excitation^ISl described by the spinon Hamil- 
tonian in Eq.(19). The only effect of the interaction U 
is to renormalizes Qf through the self-consistent equa- 
tions. The large regions of GMI indicates that lattice 
distortions of A3 > type may be helpful in the realiza- 
tion of a gapless spin liquid state in this class of materials. 
A very tiny sliver of STI is converted into a topological 
Mott insulator (TMI) above the line for which the ro- 
tors are no longer condensed. Thus, distortions of this 
type are detrimental to the realization of the TMI phase 
and suggest this phase may not be stable against lattice 
distortions. 

Next we consider the case of A3 < 0. The phase dia- 
gram is shown in the lower panel of Figj2] One distinc- 
tive feature compared to the case of positive distortion is 
that the critical Uc for the Mott transition (dashed line) 
grows as the spin-orbit coupling strength is increased. A 
second feature is that a variety of phases appear from 
the interplay of spin-orbit coupling and correlation ef- 
fects. At zero interaction, while small and intermediate 
values of spin-orbit coupling favor the metallic phase, a 
gap is opened for A ~ 3. It — 3.3t. According to the Z2 
classification'^® this insulating phase is a weak topological 
insulator (WTI) with Z2 invariant (0;010), and the small 
window survives and persists up to intermediate inter- 
actions. To the best of our knowledge, this is the first 
identification of a weak topological insulator in an inter- 
acting model. Gapless modes along certain classes of de- 
fects may be a way to identify this state in experiment.!^ 

Another interesting feature of the phase diagram for 
A3 < is that the correlation effect can drive the metal- 
lic and weak topological insulator phases into a strong 
topological insulating phase. All realized phases in the 
condensed rotor phase (below the dashed line) have car- 
riers with both spin and charge. When the correlation 
is strong enough to strip the charge degree of freedom 
(above the dashed line), those phases will turn into the 
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FIG. 2. Phase diagram of the j^ff = 1/2-band model cor- 
responding to = 5 with positive A3 = 2t (upper panel) 
and negative trigonal distortion A3 = —2t (lower panel). The 
dashed line separates the rotor condensed phases (below) from 
the uncondensed phases (above). We set t = 1, and the phases 
labeled are as follows: Strong topological insulator (STI), 
Weak topological insulator (WTI), Gappless Mott insulator 
(GMI), Topological Mott insulator (TMI), Weak topological 
Mott insulator (WTMI) and Metallic phases. 



corresponding phases with only spin degrees of freedom. 
As before, the metaUic and strong topological insulator 
phases are transformed into the GMI and TMI phases, 
respectively. Moreover, the weak topological insulator 
phases realize a novel the weak topological Mott insu- 
lator (WTMI) with increased interactions. These latter 
phases are absent on the undistorted lattic^^^ and to the 
best of our knowledge, is the first time the WTMI phase 
has been identified in a calculation. It will have gap- 
less thermal transport along the same class of defects 
that would have gapless charge (and thermal) transport 
in the WTI phase. We note that there is an acciden- 
tal gap closing in the TMI phase where two TMI phases 
are separated by a boundary. However, since the gap 
closing occurs at an even number of Dirac-like nodes, the 
topological properties remain unchanged through the gap 
closing points. 



V. j = 3/2-BAND MODEL: PHASE DIAGRAM 
OF UNDISTORTED LATTICE 

To date, the search for time-reversal invariant insula- 
tors in transition metal oxides has primarily focused on 
the j — 1/2 manifold because of its obvious connection 
to the s = 1/2 manifold heavily studied in the theoretical 
literature thus far.^^ 33 60 \vhile the j = 1/2 mani- 
fold is relevant for 1/2-filling {rid = 5), the j — 3/2 mani- 
fold is relevant for Ud — 2 which occurs in Cd2Re207 and 
Ud — 3 which occurs in Cd20s207.^^ In this section we 
investigate whether topological phases are still possible 
in the j — 3/2 manifold for some range of A and study 
the phase diagram in the presence interactions, as was 
done for j = 1/2 in Ref.jT3]. In the next section, we will 
consider the effects of distortion on the j = 3/2 phase 
diagram. 

In our calculations, we find that the non-interacting 
model with 5(P remains metallic for all physical spin- 
orbit coupling, though we can open a gap by distortion. 
A direct evaluation of the Z2 invariant shows that the 
distortion-induced insulating phase is a trivial insulator. 
So we will focus on the case with a 5cfi electr on co n- 
flguration, with Cd20s207 one possible examplc^^^^ In 
the non-interacting limit and for small values of spin- 
orbit coupling, we find a metallic phase. However, for 
A ~ 2.5t a gap opens and a STI appears. The STI phase 
is characterized by the Z2 indices (1; 000) and survives to 
moderate interaction. The corresponding phase diagram 
is shown in Fig. [3j The metallic phase is still present 
in a large portion of the phase diagram if the spin-orbit 
coupling is not too strong. It results in part from the 
fact that corresponding band structure coming from the 
j — 3/2 manifold has a larger band width than the upper 
j = 1/2 manifold and therefore has a weaker correlation 
effect for the same value of U. 

The phase boundary that separates the STI phase from 
the weakly interacting metallic phase can be determined 
analytically.'^ One notes that in the noninteracting limit 
a gap is opened at A w 2.5t. Therefore, the A and U 
that respect A ~ 2.5(5/ (C/) form the critical line between 
the STI and the metal. As in the case of the j = 1/2 
manifold, the Mott transition between the metallic (STI) 
phase and the GMI (TMI) phase is characterized by a 
vanishing amplitude of the rotor condensate at some crit- 
ical value Uc{X) (given by the dashed line). The GMI- 
metal transition and its extension to finite temperature, 
which appears to be a possible Mott transition between a 
spin-liquid insulator and metal, has already been studied 
with a possible co nnectio n to the experimental observa- 
tions in Na4lr308.'2Sl2£Ml 

The TMI phase is described by a gapped bulk spectrum 
of the spinon Hamiltonian, and is in the same topological 
phase appearing in the j — 1/2 manifold. It supports 
gapless surface states of charge-neutral spinous. The 
same spin-charge separation also occurs in 2D cases,'^^'^ 
where a quantum spin Hall state turns into an exotic 
quantum spin Hall effect at an intermediate regime of 
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FIG. 3. Phase diagram of j — 3/2-band model on the undis- 
torted lattice for = 3. The abbreviations used are the same 
as those in Fig[2] The dashed line separates the rotor con- 
densed phase (below) from rotor uncondensed phase (above). 
All energies are expressed in units of t, as before. Compared 
to the corresponding phase diagram for j — 1/2 with nd=5, 
the STI and TMI occupy a much smaller portion of the phase 
diagram.'^ The phase diagram for j — 2/3 with Ud = 2 has 
no topologically non-trivial phases within our model, even in 
the presence of distortion. 



become better defined (the self-energy scales as the en- 
ergy itself, up to logarithmic corrections).!^ Thus, the 
lowest order calculation in the U(l) gauge fluctuations 
suggests they are marginal; a more careful scaling anal- 
ysis suggest they are actually marginally irrelevant. 

Unlike the STI phase, the surface states of the TMI 
phase can not be characterized by electrical transport 
measurements due to the charge neutrality of the spinous 
on the surface of the TMI. Moreover, because this neu- 
trality, there are no Friedel oscillations around a charged 
impurity on the surface. However, spinon surface states 
of the TMI can be detected in thermal measurements, 
and by the way in which they modify the RKKY inter- 
action between magnetic impurities at the surface.l^ In 
the GMI phase, on the other hand, the bulk specific heat 
behaves^ as C ~ Tln(l/T), while in the metallic state 
it behaves as C ^ T. 



VI. j = 3/2-BAND MODEL: EFFECT OF 
DISTORTIONS 



on-site Hubbard interaction. However, the 2D nature of 
the phase suffers from an instability due to fluctuations 
of the gauge field.l^^ 

An analysis of symmetries reveals that the Hubbard 
model in Eq.([7| has a,U{l) gauge symmetry as the slave- 
rotor representation of the physical electron in Eq.(lO) 



is invariant under the following gauge transformation: 



e'^f'-fia and 



9i — (pi- In the insulating ex- 



otic state, the rotors can be integrated out since the 
charges are gapped. The resulting theory is a compact 
U{1) gauge theory coupled to the spinous.*^ The later 
theory is not stable against the fluctuations of the gauge 
field as it is a confining compact theory in 20.12^! Such con- 
finement renders the states unstable in 2D as the gauge 
fields confine the free spinon-like excitations, effectively 
removing them as legitimate low-energy excitations.'SSHl2l 
While it is believed that the extension of the spin index 
to flavors renders it deconfining for sufficiently large 
A'', the value of the critical N is not known.^il In order 
to stabilize the edge modes, the gauge fluctuations must 
be screened by other gapless degrees of freedom. In the 
2D case, this can be done by use of a bilayer structure 
in which the "second" layer contains the necessary gap- 
less degrees of freedom. In spite of the shortcomings 
of the slave- rotor mean-field theory in 2D, we note that 
recent quantum Monte Carlo calculations on the Kane- 
Mele-Hubbard model show a similar phenomenology in 
some respects at intermediate interaction strength.ESHZll 
Returning to 3D, we note that in the TMI phase the 
gapless spinon surface states are coupled to the bulk 3D 
gauge fields.!^ Thus, the low energy theory of the TMI 
phase is given by the spinon surface states coupled to the 
3D gauge fields. This theory is believed to be stable,E5Ias 
the gauge propagator is suppressed so that the spinous 



Having obtained the phase diagram of the undistorted 
j = 3/2 model in Figjsj we now study the effect of the lo- 
cal distortion of octahedra introduced in SecHIl We focus 
on two kinds of distortion with different symmetries: (1) 
a trigonal distortion of oxygen octahedra that preserves 
the C3 symmetry, and (2) a compression and elongation 
of the oxygen octahedra that preserves the C4 symmetry. 
We use parameter A3 to describe the C3 distortion and 
A4 to describe the C4 distortions. The relevant Hamil- 
tonians are given in (H and ^. 

We restrict our attention to the weak and interme- 
diate interaction limit, so we neglect possible magnetic 
phases that could become favorable in the strong corre- 
lation limit. In that limit geometrical distortions could 
alter the isotropic antiferromagnetic superexchange and 
the combined effects of spin-orbit coupling and distortion 
can give rise to an anisotropic pseudo-spin Heisenberg 
model for some perovskites.^^ 



A. Trigonal Distortion of Oxygen Octahedra 

We first consider the trigonal distortion in Eq.([2) on 
the j — 3/2, 11(1 — 3 manifold. The resulting phase 
diagrams for fixed A3 and fixed U are shown in Fig|4] 
We find the general structure of the phase diagram de- 
pends on A3 in a complicated way. As one example, 
in the upper panel of Fig|4] we show the case of fixed 
trigonal distortion A3 = 2t. In that case, we find that 
the non-interacting model is dominated by the metal- 
lic phase, which is analogous to the trigonal distortion 
driven metallic phase on the j = 1/2-band model, as 
discussed in Sec|IV| and also in Ref . [60] . The gap opens 
at A « 1.8t and closes at A « 2At. Therefore, in the 
presence of trigonal distortion only a small window of A 
admits the STI phase in the non-interacting limit. How- 
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FIG. 4. Phase diagram of the j = 3/2 band model with 
Hd — 3, including the trigonal distortion of the octahedra. 
The labeling of the phases is the same as that used in Fig[2] 
In the upper panel A3 = 2t. In the lower panel the interac- 
tion is fixed a,t U — 41 and the strength A3 of the trigonal 
distortion is varied, illustrating possible phases that may arise 
upon the application of pressure to a real system. All ener- 
gies are expressed in units of t. In both phase diagrams the 
dashed line separates the rotor uncondensed phase (above) 
from the condensed phase (below). We note the "pocket" of 
STI around A ~ 1, A3 ~ 1.5 in the lower figure has a numeri- 
cally difficult to determine boundary with the metallic phase; 
we have present our best assessment. 



erate states which in turn contribute to the formation 
of bands. However, upon the inclusion of distortion this 
degenerate manifold splits into two Kramers pairs sepa- 
rated by amount of energy related to the strength of the 
distortion, i.e. A3. The corresponding bands will also 
be separated by the same energy scale. Thus with dis- 
tortion, we are dealing with a half-filled band, with an 
effective bandwidth is reduction. So, a smaller Hubbard 
interaction is needed for the Mott transition.'iSl 

Second, distortion stabilizes the TMI phase by extend- 
ing its region of the phase diagram in comparison with 
the small region seen around [/ w 9 in the undistortcd 
lattice. (See Figjsj) We note that the GMI phase is 
found at both small and large spin-orbit coupling in the 
presence of a trigonal distortion. 

The lower panel in Fig|4] explicitly shows the effect of 
distortion at fixed interaction U — At, which is relevant 
to the application of pressure, for example. At small dis- 
tortions A3 ^ t, most of the phase diagram is dominated 
by metallic and STI phases. One can think of distor- 
tion as a driving parameter that transfers system from 
the rotor condensed phase (below the dashed line) into 
uncondensed phase (above the dashed line). Although 
the actual form of the geometrical distortion could be 
more complicated than the one we considered here, the 
result is appealing as this minimal distortion can drive 
the system across a variety of phases. Starting from the 
STI phase at zero distortion, the ground state of the sys- 
tem can exhibit a metallic behavior or perhaps transits 
to GMI and TMI phases with increased distortion. We 
hope this observation will help motivate new classes of 
experiments searching for exotic quantum phases in cor- 
related materials with strong spin-orbit coupling. 



ever, this small window grows with small but increas- 
ing interaction strength which helps to stabilize the STI 
phase. The metallic phase, however, remains dominant 
for A > 2.At, even in the presence of interaction. We note 
that the boundary separating STI phase from the metal- 
lic phase around A ~ 2.4t is not exactly but very close to 
a straight line. 

From the phase diagrams shown in Fig|4| it is clear 
that trigonal distortion has two remarkable effects. First, 
compared with the undistorted phase diagram in FigjS) it 
is evident that the distortion drives the system across the 
Mott transition (indicated by the dashed line) at rather 
smaller critical values of interaction Uc{X). This finding 
suggests that the distortion may help stabilize the TMI 
phase in a physically realistic range of interactions be- 
fore the system undergoes a transition to a magnetically 
ordered phase at strong interaction. The considerable de- 
crease of the critical Mott transition point can be traced 
back to the effect of distortion on the j = 3/2 manifold. 
Without distortion this manifold represents four degen- 



Compression and Elongation of the Oxygen 
Octahedra 



In this subsection we study the effects of the second 
type of distortion, Eq.([3]), which describes a tetragonal 
distortion of the octahedron along one of its axes. This 
distortion preserves the C4 rotation of an octahedron 
about the elongated axis, say the z-axis in Fig(l}). At 
zero spin-orbit coupling, the degeneracy of the t2g mani- 
fold will be lifted by this distortion. Compression of the 
octahedron, A4 > 0, lowers the energy of the dxy or- 
bital (which is at zero energy by our convention) below 
that of the doubly degenerate dyz and dzx orbitals, with 
energy A4. An elongation (expansion), A4 < 0, of an 
octahedron lowers the energy of the dyz and dzx orbitals 
relative to dxy This rearranging orbitals strongly affects 
the magnetic properties of the double perovskites in the 
strong interaction limit.'^ 

When spin-orbit coupling is present, the levels split in 
a more complicated way. Similar to Eq.(29) 



for the e'g 

manifold, the spin-orbit coupling results in the following 
effective Hamiltonian for a proper linear combination of 
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\yz) and \zx) states^ 



{H.. 



SO/ yz^zx 



(30) 



the doublet {~^i\v^) + 



and a is the usual Pauli ma- 



where r acts within 

Azx)), ■^{\yz) - i\zx))}, 
trix of real spin. Note that the spin-orbit coupling acts 
like a Zeeman coupling so that the effective magnetic 
field has opposite direction in different states of the dou- 
blet, therefore, the time-reversal symmetry is preserved. 
[Evident as well from the Hamiltonian ([t]).] 

The effect of distortion on the spin-orbit basis can also 
be treated in the same way. In particular, we consider 
its effect on the quadruplet j = 3/2 manifold. (See Ap- 
pendix [b]) With distortion the following states are ob- 
tained: 



=--^(|yz ^)+i\zx^)), 



with energy £1^2 = |(2A4 - 

IV'3) -C^[-/(A,A4)|2;z ;> 
1^4) =C[/(A,A4)|2;z t> 



A) and 

~if{\,^^)\zx t> 



(31) 



\xy i)] , 
(32) 



with energy £3,4 ^ \{'2.l^i + A - yiAjTlA^ATQA^ ) . 
Here C is a normalization constant depending on /. Note 



that in Eq. ( 32 1 / is a function of its arguments with 
/ — > 1 as A4 — 0. In the limit of vanishing distortion, 
the above sates IV'1,2) and IV'3,4) reduce to the four sates 
||,±|) and ||,±i) of the quadruplet j = 3/2 mani- 
fold, respectively. Note the |V-'i,2) keeps the character of 

= if states even for A4 ^ 0. With electron occupa- 
tion Ud ~ 2> and for A4 > (A4 < 0), the states \ipi^2) 
(l''/'3,4)) form a half filled band, and we will see that the 
C4 distortions strongly affect the phase diagram found in 
Fig|3]for the undistorted lattice. 

Figjs] depicts the phase diagram of the j = 3/2 model 
with compression (elongation) distortion in upper (lower) 
panel. While both compression and elongation of octa- 
hedra possess almost the same critical (dashed) line for 
the Mott transition, the topological phases occupy rather 
different regions. For example, the case of compression 
does not support a STI phase, while elongation does. 
For strong compressional distortion A4 > 0, the bands 
are mainly comprised of the states IV'1,2), which are not 
spin-orbital entangled. The states IV'3,4): on the other 
hand, are spin-orbital entangled states. Even for elon- 
gation, most partions of the condensed phase (below the 
dashed line) is dominated by the metallic phase even in 
the presence of interaction. 

For both signs of the C4 distortions, strong interac- 
tions open up a gap in the bulk spinon spectrum turning 
the GMI phase into the TMI phase as shown Figj5] For 
A4 > 0, a weaker distortion may extend the boundary 
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FIG. 5. Phase diagram of the j = 3/2-band model with C4 
tetragonal distortion. The upper and lower panel correspond 
to compression and elongation distortion of octahedra, re- 
spectively. We set A4 = 2t for compression and A4 — ~2t 
for elongation. All energies are expressed in units of t. In 
both phase diagrams the dashed hne separates the rotor un- 
condensed phase (above) from the condensed phase (below). 
Note that for these values of distortion, the topological phases 
are interaction-driven, i.e. they do not extend down to the 
non-interacting hmit as they do in Figs. 2-4. 



of the TMI phase towards weaker interactions, and per- 
haps open a region with a STI phase. This is because at 
weaker distortion there would be a considerable contri- 
bution from unentangled states in the formation of the 
bands. This is clearly seen in the lower panel of Figj5] 
where a finite region with a STI phase is established. If 
the spin-orbit coupling is kept fixed, at very strong distor- 
tion the coefficient / in Eq.(32) tends to zero, and there- 

This may partly 



fore the states become unentangled. 
explain why at small interaction the metallic phase is 
dominant. However, the lower panel of Fig js] reveals that 
interactions can drive the formation of strong topolog- 
ical insulators even when the STI phase is not present 
in the non-interacting limit, as it is in most of the cases 
considered previously. A two dimensional analogue of 
this problem has been studied elsewhere,I^SHIII where it is 
shown that the interaction-driven insulating phases with 
nontrivial topology can be found on a variety of different 
lattices. 

While the distortion favors the metallic phase in the 
weak interaction limit and the GMI phase in the strong 
limit (above the Mott transition) compared to the non- 
distorted case (Figj3|, it is possible that strong dis- 
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order can transform the metallic phase and GMT into 
topological phases, the so called topological Anderson 
insulator PSHSll particular, the effect of disorder on the 
GMI phase would be an interesting problem. 



VII. SUMMARY AND CONCLUSIONS 

In this work we investigated the phase diagram of some 
transition metal oxides with 5d orbitals on the pyrochlore 
lattice. We focused on the interplay between electron cor- 
relation, spin-orbit coupling, and distortion. Our main 
results are summarized in the phase diagrams presented 
in Figs. [2][5] obtained within the slave-rotor mean- field 
theory. 

Examples of pyrochlore transition metal oxides include 
A2lr207 (A=Y, Pr, Eu or other rear-earth elements), 
Cd20s207, and Cd2Re207, in which different transition 
ions favors either j = 1/2 or j = 3/2 manifolds to be 
partially occupied .'^^'^ A central feature of our work was 
to considered distortion of local oxygen octahedra sur- 
rounded the transition ion in the presence of interac- 
tions. Such distortions are inherent to the systems we 
studied.EHESl 

We first studied the effect of trigonal distortion on the 
j = 1/2 phase diagram already obtained in Ref.|15j. In 
the noninteracting limit, we found the distortion destabi- 
lized the STI phase and turned it into a metallic phase.l^ 
(See Figj2j) However, a A3 < can also help stabilize a 
weak topological insulating (WTI) phase which becomes 
a weak topological Mott insulator (WTMI) above the 
Mott transition line. To the best of our knowledge, these 
features have not been obtained in previous interacting 
models before. 

We also extended the study of interacting topologi- 
cal insulators in transition metal oxides to include the 
case where the j = 3/2 manifold is partially filled. (See 
Figjs]) We found that strong spin-orbit coupling opens a 
gap in the non-interacting spectrum and the STI appears, 
along with a metallic phase at small spin-orbit coupling. 
These phases persist in the presence rather large inter- 
actions interaction due to the large band width of the 
non-interacting model, and eventually at large enough 
interactions the Mott phases appear. Most portions of 
the Mott phase are identified as a gapless Mott insulator 
(GMI). However, at some intermediate regime of spin- 
orbit coupling, 1 ^ A ^ 2, the TMI phase is obtained. 

Trigonal distortion extends the TMI phase to a wider 
range of interaction and spin-orbit coupling. (Fig|4]) 
Trigonal distortion also decreases the critical interaction 
for the Mott transition. Moreover, we showed that the 
distortion can serve as a tuning parameter in which the 
transition between a variety of phases could occur by dis- 
torting the lattice, though a more realistic form of dis- 
tortion could have a more complicated evolution of the 
phases. We also examined the effect of tetragonal distor- 
tion of octahedra caused by an elongation or compression 
of an octahedra along one of its axis. (Figjs)) The STI 



phase is found to be very delicate with respect to this 
type of distortion, and most of the phase diagram is occu- 
pied by either the metallic phase the the GMI. For strong 
enough distortion, however, interaction can restore both 
STI and TMI phases. The restoration of these phases is 
an example of "interaction-induced" topological phases 
as these phases do not persist down to zero interaction. 

One might wonder to what extent the slave-rotor 
mean-field results should be trusted. Is there an alter- 
native method that can be used to obtain a TMI phase, 
for example? As we mentioned earlier, quantum Monte 
Carlo methods applied to the 2-d Hubbard model on the 
honeycomb lattice seem to suggestl^^HZH that there is an 
intermediate, gapped phase that lives over a similar re- 
gion of the phase diagram that the slave-rotor method 
predicts a 2-d TMI (recall that the slave-rotor method 
is not expected to be reliable in 2-d). This may suggest 
that there is indeed a state with fully gapped bulk exci- 
tations, but with gapless spin excitations on the bound- 
ary. However, it may be that the slave-rotor method 
fails to correctly capture the collective nature of the 
"true" low-energy spin-excitations by forcing them into 
a single-particle mean-field formalism. One may also ask 
about the reliability of the slave-rotor method more gen- 
erally. The original work of Florens and Georgea^SED 
on Hubbard models shows a favorable comparison with 
dynamical mean-field theory and Gutzwiller projection 
for quantities like the quasi-particle weight and effec- 
tive mass below and just above the Mott transition. Fi- 
nally, the references contained in the work of Pesin and 
BalentJi^ provide further support for the reliability of the 
slave-rotor method when compared with path-integral 
renormalization group calculations and variational clus- 
ter methods on frustrated lattices. Taken together, it 
seems the method works reasonably well in situations 
where interactions are not too strong and no magnetic 
order is expected. Nevertheless, a more careful study of 
the possibility of a TMI phase within a more sophisti- 
cated class of calculations remains highly desirable and 
we hope this work will help to inspire such studies. 

Regarding the physics of pyrochlore oxides, a number 
of interesting directions for future study remain. For ex- 
ample, it would be highly desirable to have a better un- 
derstanding of the specific form of lattice distortions that 
occur in nature and what their influence is in terms of 
candidate topological phases. It would also be interest- 
ing to obtain a better understanding of disorder on the 
interplay of correlations, spin-orbit coupling, and lattice 
distortions. Finally, we note that even more exotic possi- 
bilities exist for novel phases when certain conditions are 
j^g^pBHHSI improved understanding of how likely the 
conditions for these "fractional" phases with non-trivial 
ground state degeneracy are to be met in real materials 
would be welcome. 

As we restricted ourselves to interactions that were 
not too strong, we did not invoke the possible mag- 
netic phases that could be more favorable at very strong 
Coulomb interaction. The magnetic phase is interest- 
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ing in its own right as the pyrochlore lattice has a ge- 
ometrically frustrated structure. The latter property 
along with the proximity to the metallic phases can 
enhance the quantum fluctuations. Hence, even the 
nearest-neighbor antifcrromagnetic interaction may sta- 
bilize a spin-liquid phase on the pyrochlore lattice.'^SHIIl 
Besides the antifcrromagnetic interaction, j = 1/2 mag- 
netic models that include some additional interactions 
such as Dyzaloshinsky-Moriya and other anisotropic in- 
teractions can help sustain orde ring on the pyrochlore 
lattice at low temperatures.'2212ll However, the situation 
is more complicated for the j = 3/2 model: because of 
orbitally-dependent exchange, biquadratic (forth order 
in spin operators) and triquadratic (sixth order in spin 
operators) interactions arise.*^ These new interactions 
give rise to some exotic phases in double perovskites,!^ 
and tetragonal distortion of octahedra can result in 
an anisotropic pseudo-spin antifcrromagnetic exchange 
Heisenberg model.l^ Such models can be developed for 
our model with distortion, too. Indeed, what magnetic 
phases become favorable and how they are related to the 
topological phases we addressed here are interesting open 
problems. 



Appendix A: Transformation to orbital angular 
momentum states 

In terms of the t2g states, the effective I = 1 angular 
momentum states are given by 



1^ 



' V2 
\l, = 0) = \xy), 



{\yz) + i\zx)), 



\iz = -1) = " *l^2;)), 

so that in terms of the orbitals the effective Iz angular 
momentum is given by 




in the basis {dxy,dyz,dzx)- This immediately gives 
+ in terms of the state occupations. 



Appendix B: Spin-orbit coupled states 



The transformation between the spin s = 1/2 in the 
effective / = 1, t2g orbital basis {dxy, dyz, dzx) and the 
basis of \j,jz} for j — 1/2, 3/2 is given by 
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